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ABSTRACT 

We describe radio observations at 244 and 610 MHz of a sample of 20 luminous and 
ultra-luminous IRAS galaxies. These are a sub-set of a sample of 31 objects that 
have well-measured radio spectra up to at least 23 GHz. The radio spectra of these 
objects below 1.4 GHz show a great variety of forms and are rarely a simple power- 
law extrapolation of the synchrotron spectra at higher frequencies. Most objects of 
this class have spectral turn-overs or bends in their radio spectra. We interpret these 
spectra in terms of free-free absorption in the starburst environment. 

Several objects show radio spectra with two components having free-free turn- 
overs at different frequencies (including Arp 220 and Arp 299), indicating that syn- 
chrotron emission originates from regions with very different emission measures. In 
these sources, using a simple model for the supernova rate, we estimate the time for 
which synchrotron emission is subject to strong free- free absorption by ionized gas, and 
compare this to expected HII region lifetimes. We find that the ionized gas lifetimes 
are an order of magnitude larger than plausible lifetimes for individual HII regions. 
We discuss the implications of this result and argue that those sources which have a 
significant radio component with strong free-free absorption are those in which the 
star formation rate is still increasing with time. 

We note that if ionization losses modify the intrinsic synchrotron spectrum so that 
it steepens toward higher frequencies, the often observed deficit in fluxes higher than 
^10 GHz would be much reduced. 
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1 INTRODUCTION 

The radio emission from star forming galaxies is the result of 
two emission processes, non-thermal synchrotron, and ther- 
mal Bremsstrahlung, or free-free emission. Synchrotron is by 
far the dominant component in the 1-10 GHz range and is 
characterized by a power-law emission spectrum, /„ oc v~ a 
where a ~ 0.8. Free-free emission has an almost flat spec- 
trum and the emission falls with frequency only as i/ -0 ' 1 . 
Therefore, towards higher frequencies free-free emission be- 
comes increasingly important and radio spectra are expected 
to show a flattening beyond ~ 10 GHz. 

In LIRGs and ULIRGs, although such a flattening is 
sometimes seen, the majority of radio spectra, even when 
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measured to 23 GHz, show no such flattening. In fact, there 
is a tendency for the spectra to steepen (Clemens et al. 2008). 
Despite this, Clemens et al. (2008) showed that the spectral 
index between 8.4 and 23 GHz is flatter for sources with 
high values of the far-infrared-radio flux ratio, q. As these 
sources, with an excess of infrared relative to radio emis- 
sion, are expected to be younger, they would be expected to 
have proportionately more free-free gas, and thus shallower 
spectral indices at high frequencies. 

As we will see later, the radio spectra of LIRGs and 
ULIRGs also hold surprises towards low frequencies, with 
spectra showing a variety of bends and turn-overs. From a 
theoretical standpoint there is more than one way in which 
the radio spectrum may flatten towards low frequencies. 
There are those processes that absorb radio photons, free- 
free absorption and synchrotron self-absorption and those 
that cause energy losses of the relativistic electron popula- 
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Table 1. Total 244 and 610 MHz fluxes. Distances are as reported 
in Condon et al. (1991). 
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tion directly. Of these, ionization losses have been considered 
by Thompson et al. (2006) as a cause for the flattening of 
starburst radio spectra towards low frequencies. However, as 
illustrated by these authors, the flattening is rather gradual. 
Our data show evidence of rapid changes in spectral index 
and also inversions of the spectra, much more consistent with 
an absorption process with a strong energy dependence. 

Both free-free absorption and synchrotron self- 
absorption could rapidly absorb the power-law synchrotron 
spectrum toward low frequencies. However, synchrotron self- 
absorption requires extremely high brightness temperatures 
in order to be significant, and these are just not attained 
(with the possible exception of UGC 8058) in our sources 
(Condon et al. 1991). We therefore consider free- free absorp- 
tion to be the defining mechanism for the spectral shapes of 
this sample, especially at low radio frequencies. 

Here we describe new observations of a sample of LIRGs 
and ULlRGs at 244 and 610 MHz made with the GMRT, 
which we combine with archival data and our previous 
23 GHz data to produce well-sampled radio spectra across 
this wide frequency range. We consider the structure of the 
starbursts in these objects that is implied when their radio 
spectra are described in terms of the free-free absorption of 
a synchrotron power-law. 



2 SAMPLE SELECTION 

Our sample is based on that of Condon et al. (1991) who 
made 8.4 GHz observations of the 40 ULIRGs brighter than 
5.25 Jy at 60 fim. In Clemens et al. (2008) we described 
new observations at 23 GHz and the collection of archival 
data for a sub-sample of 31 of the 40 described by Condon 
et al. (1991). Some objects were excluded for observational 
reasons, as described in Clemens et al. (2008). Only 20 of 
these 31 objects were observed at the GMRT at 244 and 
610 GHz simply due to observing constraints. 



3 OBSERVATIONS 

Observations at 244 and 610 MHz were made between 2008 
January 13 and 16 at the Giant Metre-wave Radio Tele- 
scope, (GMRT). The two frequencies were observed simul- 
taneously using the 'dual' mode, with one polarization for 
each frequency. The total bandwidth of 8 (16) MHz was split 
into 64 (128) channels at 244 (610) MHz. Total integration 
times for the sources were between 30 and 120 minutes with 
scans on-source being interleaved between scans of phase 
calibrators approximately every 30 minutes. Both 3C 147 
and 3C 286 were used as flux calibrators. 

Data reduction was carried out using the NRAO As- 
tronomical Image Processing System AIPS. After initial, 
channel-based, editing to remove data affected by interfer- 
ence, bandpass calibration was applied using one of the flux 
calibrators. The data were calibrated in the standard way 
by forming a "channel 0" data set from the average of 10 
channels near the centre of the band. The complex gains 
from this process were then copied to the line data before 
reducing the number of channels by a factor of 10 for the 
purpose of map making. 

Maps were made by Fourier transformation of the data 
in each channel to produce 'dirty' images and beams (aips 
task "imagr"). The images and beams were then combined 
to produce a single continuum image and beam for each 
source, thus avoiding bandwidth smearing. The resulting 
maps were deconvolved using the clean algorithm and the 
combined beam. Where the source was sufficiently bright, 
self-calibration was further applied (varying only the an- 
tenna phases) to improve the signal-to-noise. The final maps 
typically had resolutions of 13" and 5" and rms noise levels 
of 5 mjy and 0.5 mjy at 244 and 610 MHz respectively. 

Flux densities were measured by direct integration on 
the maps via the MPS task tvstat. These integrated fluxes 
(reported in Table 1) have errors which include both the rms 
noise and calibration errors. 



4 RESULTS 

4.1 Low frequency radio spectra 

The radio spectra shown in figure. 1 show various forms 
at low frequencies. These include an almost straight power- 
law from 22 GHz down to 244 MHz (e.g. IC 1623), a sin- 
gle turn-over at low frequencies (NGC 1614), evidence of 
two turn-overs (e.g. IR 03359+1523) and a gradual flatten- 
ing towards low frequencies (e.g. NGC 2623). In addition, 
IR 08572+3915 shows an almost flat spectrum. This may rise 
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Figure 1. Model fits to the radio data. The radio fluxes used, other than those given in table 1, are collected in Clemens et al. (2008). 
The turn-over frequencies are indicated with vertical lines. Poor quality fits in which the adopted model seems insufficient are shown as 
dot-dashed lines. 



at frequencies below 1 GHz but the data point at 610 GHz 
has a large error and is marginally consistent with a spec- 
trum that remains flat. 

The flattening of the radio spectra is due to absorp- 
tion of the power-law synchrotron emission, S v oc v~ a , 
where v is frequency. In all but one case this must be free- 
free absorption by ionized gas in the sources because the 
brightness temperatures are not high enough for synchrotron 
self-absorption to be important. Only UGC 8058 1 has an 
8.4 GHz surface brightness sufficiently high that synchrotron 
self-absorption may be important (Condon et al. 1991). 

The shape of the low-frequency spectrum depends, al- 
most entirely, on the quantity and distribution of ionized 

1 Although this source is known to contain compact radio emis- 
sion from an AGN which shows variability on timescales of years, 
Ulvestad et al. (1999) show that this variability is likely confined 
to frequencies above 5 GHz. In our spectrum the data points at 
1.4, 4.8, 8.4, 15 and 22.5 GHz were, in any case, all made on 
the same date and a turnover is already seen between 1.4 and 
4.8 GHz. Sec section 4.4. 



gas with respect to the source of synchrotron emission. The 
free-free optical depth is, t„ oc n 2 Lv~ 21 , where n e is the 
ionized particle density, v is the frequency, and L is the path 
length along the line-of-sight. n 2 L is the emission measure. 
Radio spectra will show a turn-over where t„ = 1. If the 
frequency of this turn-over is vt then the optical depth is 
given by, r„ = {vjvtY 2 ' 1 . 

The radio spectra of normal, star-forming galaxies have 
been reviewed in detail by Condon (1992). The spectra are 
the result of emission from a medium containing both syn- 
chrotron and free-free emitting gas. Based on observational 
results, Condon assumed that the thermal emission con- 
tributes 10 times less than synchrotron emission at 1 GHz. 
Based on equation 6 of Condon, the radio surface brightness 
therefore varies with frequency according to, 

where S is flux density, f2 is the area subtended by the 
source and T e is the electron temperature. H is the ther- 
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Figure 1. Continued 



mal fraction. For synchrotron emission that contributes 10 
times that of thermal emission at 1 GHz, as assumed by 
Condon (1992), H = 0.1. Fig. 4 of Condon shows the family 
of curves that result from varying the turn-over frequency, 
the frequency at which the optical depth t v = 1. In fig. 2 
we illustrate the effect of varying H for a fixed turn-over 
frequency. 

Because we do not have measurements of the the source 
size, O, at every measured frequency equation 1 is only ap- 
plicable under the assumption that the size of each source is 
not a function of frequency. Although cosmic ray diffusion is 
known to make this untrue for some nearby objects, the very 
short synchrotron lifetimes expected in the present sample 
should ensure very similar source sizes at all frequencies un- 
der study here. 

For normal galaxies the spectral index is close to —a 
at intermediate frequencies (5 GHz) and flattens towards 
higher frequencies due to the increasing contribution of free- 
free emission. Towards low frequencies the spectra turn-over 
due to the effects of free-free absorption. The frequency at 
which this turn-over occurs depends on the emission measure 
of the ionized gas. As pointed out by Condon, real sources, 
such as M 82, may be expected to show a flattening of the 



radio spectrum towards low frequencies rather than a turn- 
over due to source inhomogeneity. 

We use equation 1 as a starting point for the analysis 
of our radio spectra of luminous IRAS galaxies. 

Each turn-over frequency in our spectra, and some have 
more than one, is the result of absorption by ionized gas 
with a different emission measure. Therefore, by fitting one 
or more components described by equation 1 we can de- 
termine L for our sample galaxies via the standard ex- 
pression for the free-free absorption coefficient (equation 1 
in Condon 1992). This step requires an assumption for the 
electron temperature, which we take to be 10 4 K. In fitting 
equation 1 to the data we consider the following cases. 

(i) The spectrum is straight to the lowest measured fre- 
quency (e.g. IC 5298)- no turnover can be estimated. 

(ii) The spectrum shows a single turn-over. Either only 
a single absorbed component is visible (e.g. NGC 1614) or 
there is evidence of another component whose flux continues 
to rise to the lowest measured frequency (e.g. UGC 5101). 
In this case we fit equation 1 allowing a, H and the turn- 
over frequency to vary. If data points are insufficient we fix 
q = 0.8. We allow for the presence of an unabsorbed com- 
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Figure 1. Continued 



ponent by fitting an additional power-law with the same 
power-law index. For those objects where only one compo- 
nent is seen we provide a limit to the possible contribution 
of an unabsorbed component by finding the maximum con- 
tribution of a straight power-law that is consistent with the 
lowest frequency data point. We refer to this as model 1. 

(hi) The spectrum shows 2 clear turn-overs (e.g. 
IRAS 03359+1523). We fit 2 components obeying equation 1 
but fix H = 0.1 and a — 0.8. For each component, we fit the 
turn-over frequency and the flux ratio of the components. 
If equation 1 is represented by /» then fit / = /i + x fi- 
In this case, if component has the lower break frequency 
the constant x represents the relative area covered by /2, 
ie. Q.2/0.\ because we measure integrated fluxes rather than 
surface brightnesses. We refer to this as model 2. 

Equation 1 was fit to the data via a non-linear least 
squares procedure that minimized Chi-squared between the 
data and model. We limited the number of parameters that 
were allowed to vary so as to remain fewer than the num- 
ber of data points. In table 2 we list the model used and 
associated fit parameters. In some cases, where parameter 
values were fixed, reasonable fits could not be obtained with 
a power-law index of a = 0.8 and/or H = 0.1. In these, 
alternative fixed values were used, as indicated in table 2. 

If the best fit solution passed though less than 3/4 of 
the data points, the fit was deemed to be poor. In these 
cases, shown as dashed lines in fig. 1, the derived quantities 
are uncertain and appear in parentheses in table 2. In these 
cases, a good fit would require either a more complex model, 
or more parameters to be varied within the adopted model, 
beyond those permitted by the data. 



The model fits are shown in fig. 1. We note that the val- 
ues of H, and to a lesser extent a, are often rather poorly 
constrained by the data. H tends to be poorly determined 
because it is mostly constrained by high frequency data 
points beyond 8.4 GHz. Here the fractional thermal emis- 
sion should be significant (in normal galaxies ~30% of the 
emission at 8.4 GHz is thermal). Equation 1 is a spectrum 
that flattens towards higher frequencies. However, as already 
noted in Clemens et al. (2008) the radio spectra of LIRGs 
and ULIRGs often do not show the flattening expected due 
to free-free emission and frequently show a steepening. In- 
deed, Fig. 1 has several examples where the model over- 
predicts the flux at > 10 GHz. The fitted thermal fractions, 
which are generally much lower than 0.1, are a consequence 
of the fact that only a spectral flattening is permitted by 
equation 1. This shows that there is an effect that is not 
taken into account in the simple model that we apply here 
and is discussed in section 6.1. 



The spectral index, a, on the other hand, is poorly con- 
strained by the limited number of data points when the turn- 
over frequency occurs at frequencies so high that no part of 
the spectrum is approximately straight. 



However, we are concerned mainly with identifying the 
frequency of the turn-over seen towards longer wavelengths, 
and the relative strength of the two emission components 
where a second is present. These values are only very weakly 
affected by the value of H and weakly by the value of a. 
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Table 2. Fitted model parameters. Values that are very uncertain are shown in parentheses. Derived star formation rates refer to the 
mass range 0.15-120 Mq. * parameter held fixed at this value, t UGC 8058 contains a radio loud AGN that emits above 4.8 GHz, see 
section 4.4. 
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4.2 General form of the radio spectra 

In almost all cases, our low frequency radio data have re- 
vealed the presence of a turn-over in the radio spectrum. 
Although the turn-overs are often at frequencies above 
610 MHz, such turn-overs were previously detected only at 
a single frequency. For an object with only 3 or 4 measured 
fluxes there was little justification for fitting anything other 
than a straight line to the data. 

The addition of low frequency radio data to the spec- 
tra of this sample of LIRGs and ULIRGs shows the spec- 
tra are frequently more complex than previously thought. 
Galaxies that showed evidence of turn-overs above 1.4 GHz, 
and would have been expected to show fluxes that fell with 
decreasing frequency, often show spectra that actually rise 
again (e.g. NGC 2623). Although Condon (1991) pointed out 
that free-free absorption in real sources is likely to result in 
rather broad turn-overs or just a flattening of the spectrum 
due to source inhomogeneity, we often see rather more com- 
plex spectra, indicative of at least 2 emission components. 
This was unexpected. 

Given our limited sensitivity at 244 and 610 MHz as 
well as the fact that there are rarely data below 244 MHz it 
is likely that all sources in our sample have more than one 



emission component; we either don't detect the emission or 
have not observed at sufficiently low frequencies. 

The presence of radio emission components with dif- 
ferent characteristic turn-over frequencies tells us that there 
are populations of relativistic electrons in environments with 
diverse emission measures. There is more than one physical 
scenario consistent with this result, and we will discuss this 
issue in section 5.1, but the most simple explanation may 
be that there are starbursts occurring in two or more nuclei 
(as most of these galaxies are probably merger remnants 
this would seem plausible). If either the mean ionized gas 
density or the path length through the ionized gas to the 
supernovae were different in the two nuclei, different turn- 
over frequencies would result. Reference to the 8.4 GHz, 0'.'25 
resolution maps of Condon et al. (1991) shows that although 
some sources with 2 emission components do indeed show 2 
spatially resolved nuclei (e.g. Arp 220) there are others that 
do not (e.g. NGC 2623) and also sources that show a single 
component radio spectrum but have 2 radio-bright nuclei 
(e.g. IRAS 14348-1447). Of course, higher resolution data 
may reveal additional double nuclei, and additional radio 
data at lower frequencies may reveal second emission com- 
ponents where only one is currently detected. Therefore, we 
cannot, at present, identify separate emission components 
with separate starburst nuclei. 
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Table 3. Emission measures. 



Source 


EMi 


ni 


EM2 


712 




10 6 cm~ 6 pc 


— 3 
cm J 


10 6 cm -6 pc 


— 3 
cm J 


NGC34 


<3 








IC1623 


23 


330 






CGCG436-30 


<3 








101364-1042 


11 


290 


_ 


_ 


IIIZw35 


8.2 


310 






UGC2369 


0.8 


100 


_ 


_ 


103359+1523 


140 


1400 


0.32 


68 


NGC1614 


0.24 








105189-2524 


6.2 


200 






NGC2623 


9.9 


290 






108572+3915 


56 


780 






UGC4881 


0.32 


47 






UGC5101 


3.7 


190 






110173+0828 


90 


^1200 


1.8 


>166 


110565+2448 


0.40 


35 






Arpl48 


0.20 


42 






UGC6436 


0.66 








Arp299 


62 


900 


0.38 


70 


112112+0305 


15 


310 


0.24 


40 


UGC8058 


53 


^950 


0.12 


>44 


UGC8387 


31 




0.36 




NGC5256 


0.074 


24 






UGC8696 


3.1 


130 






114348-1447 


3.1 


120 






IZwl07 


31 


460 






115250+3609 


16 


>500 






Arp220 


13 


410 


0.28 


59 


NGC6286 


0.14 








NGC7469 


13 


560 






IC5298 


<3 








MRK331 


<3 









For all objects with sufficient data we provide the 
8.4 GHz flux ratio of the 2 model components in table 2, or 
an upper limit where only a single component is detected. 



4.3 Star formation rates 

In order to estimate the star formation rates associated with 
each emission component we make use of the calibration of 
Panuzzo et al. (2003) (based on Bressan et al. 2002). This 
relates the 8.4 GHz luminosity to the star formation rate 
for stars with a Salpeter IMF and masses in the range 0.15- 
120 M via SFR = 258.4 x 10 - 23 L 8 .4 Moyr^CWHz- 1 )- 1 . 
In table 2 we give the observed 8.4 GHz luminosities and 
implied star formation rates. For those objects with 2 emis- 
sion components we show the star formation rate for each. 
For galaxies in which 2 emission components are clearly de- 
tected, both imply high star formation rates. Although a 
high star formation rate is a requirement for detection at 
frequencies below 1.4 GHz if the radio emission is derived 
from supernova explosions, the significant fraction of sources 
in which 2 such components are detected (13/23 with data 
below 1.4 GHz) means that this is not uncommon among 
such galaxies. The radio emission at low frequencies is not 
therefore from a small number of supernovae occurring in 
more quiescent regions beyond the main starburst, but is 
rather from a population of supernovae of similar size tak- 



ing place in a lower density medium (with a lower emission 
measure) . 

In those objects where only a single component is de- 
tected, the upper limits on the fractional contribution of a 
component with a lower frequency turn-over only constrains 
the star formation rate to be <10% that of the main com- 
ponent in 3 cases. The radio emission from infrared lumi- 
nous galaxies is therefore consistent with high star forma- 
tion rates in regions with 2 distinct emission measures (but 
see section 5). 

4.4 UGC 8058 

UGC 8058 (Mrk 231) is the only source in our sample with a 
radio loud AGN. Vega et al. (2008) found that the contribu- 
tion of the AGN to the bolomctric luminosity was 24%. The 
AGN is also variable, although the variability is seen only 
above 5 GHz (Ulvestad et al., 1999). Fig. 1 shows that our 
model succeeds in fitting the radio data for this source per- 
fectly well as two absorbed components. The discontinuity 
is probably due to radio emission from the AGN that is seen 
only above 4.8 GHz due to strong absorption. In this case the 
absorption may be synchrotron self-absorption or free-free 
absorption. The spectrum implies that, at the epoch of ob- 
servation of the data above 1.4 GHz (1995 November 17), the 
AGN contributed 66% to the radio emission at 8.4 GHz (see 
table 2). The radio emission at lower frequencies is therefore 
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derived almost entirely from star formation, and we can use 
this fact to estimate the star formation rate. In table 2 we 
give the star formation rate for this component alone. The 
models of Vega et al. (2008) assumed that AGN did not emit 
in the radio. Interestingly, the model for this source fitted 
well the low frequency radio data, but under-predicted the 
radio emission above 5 GHz. 

4.5 Emission measures 

Table 3 shows the emission measures derived from the fitted 
turn-over frequencies listed in table 2. Also shown are the 
mean ionized gas densities implied by calculating the phys- 
ical source sizes via the angular sizes measured at 8.4 GHz 
and distances given by Condon et al. (1991). Our observa- 
tions rarely resolve the sources so we do not have a direct 
measure of the source sizes at 244 and 610 MHz. These could 
, in principle, be larger than at 8.4 GHz due to cosmic ray 
diffusion, but the synchrotron lifetime is anyway so short 
that the effect on source size should be negligible. We can- 
not exclude the possibility, however, that the low-frequency 
radio emission comes from a more extended distribution of 
supernovae. However, such an extended component would 
also emit at 8.4 GHz and should have been detected by Con- 
don et al. as extended halo emission. It is therefore probably 
reasonable to use the source sizes as measured by Condon 
et al. (1991). 

Free-free absorption towards individual supernova rem- 
nants has been used to derive emission measures in M 82 
(Wills et al. 1997). A large fraction have turn-overs that 
are evident at 1.4 GHz, implying emission measures ^ 
10 6 cm" 6 pc. Other studies have identified regions with 
emission measures as high as 4 x 10 6 cm -6 pc (Seaquist, 
Bell & Bignell, 1985) and Clemens & Alexander (2004) de- 
rived a value of 10 7 for the centre of the luminous IRAS 
galaxy, UGC 83 87 2 , based on high resolution radio data. 
Emission measures as high as 1.8 x 10 8 cm -6 pc have been 
found for compact sources in the "Antennae" galaxies by 
Neff & Ulvestad (2000). For those objects which were ob- 
served at 244 and 610 MHz we find emission measures in 
the range 0.07 - 140 x 10 6 cm" 6 pc. 

For those sources in which 2 absorbed components are 
seen we have fixed the thermal fraction of both so that 
H = 0.1. In this case, the spectra imply that the less ab- 
sorbed component covers a much larger area than that of the 
more absorbed component, because equation 1 implies lower 
brightness temperatures for lower turn-over frequencies. Be- 
cause we measure integrated fluxes, a larger source area is 
needed to fit the data. Of course, this may not reflect the 
true source structure. The component whose turn-over lies 
at a lower frequency may instead have a much lower thermal 
fraction. The similarity between decreasing the thermal frac- 
tion by a factor of, say 10, and increasing the source size by a 
factor of 10 for a given turn-over frequency is shown in fig. 2. 
At frequencies near the turn-over, and below, the resulting 
spectra are rather similar. At higher frequencies there is a 
significant difference, as the thermal emission contributes a 
larger fraction of the total. There is little difference, then, 
for the low frequency end of the model fits between these 

2 The 2 values we derive here for UGC 8387 straddle this value. 




0.1 1.0 10.0 100.0 

frequency/GHz 

Figure 2. Comparison between varying the thermal fraction and 
the source size for a given turn-over frequency. The black lines 
are for a fixed source size and thermal fractions, H, in equation 1 
of 0.2, 0.1, 0.05, 0.02 and 0.005, from bottom to top. In this case, 
increasing flux corresponds to increasing surface brightness. The 
coloured line (top line in plot) has a thermal fraction of H = 0.1 
and has been scaled by iJ/0.005. This corresponds to increasing 
the source area by a factor of 20 at fixed thermal fraction. The 
surface brightness remains constant in this case. 



2 possibilities, although at higher frequencies we limit our 
ability to fit the data. However, as mentioned above, the 
high frequency data are poorly modeled, and so constrain- 
ing the thermal fraction at v > 10 GHz is unreliable in any 
case. 

To summarize, the relative strength of the two emission 
components either places a constraint on the relative phys- 
ical area of each (thermal fraction of both fixed) or places 
a constraint on the difference in thermal fraction between 
the two. A combination of these two extremes is also per- 
mitted, of course. In table 2 we give the value of x, which 
is the relative physical area covered by the less absorbed 
component, or, almost equivalently, the factor by which this 
component is poorer in thermal emission. For example, the 
component with a lower frequency turn-over in Arp 220 is 
either more spatially extended by a factor of 40 than that 
of the component with the higher frequency turn-over, or it 
has a thermal fraction 40 times lower. 



5 IMPLICATIONS FOR STARBURST 
EVOLUTION 

As discussed in Condon et al. (1991), with the exception of 
UGC 8058, there is little evidence that any of the galaxies in 
our sample contain a radio loud AGN. Vega et al. (2008), in 
fact, demonstrate that star formation alone is sufficient to 
explain the radio emission of every source except UGC 8058. 
In what follows, we assume that the radio emission is derived 
entirely from processes linked to star formation. 

In this section we focus particularly on the presence of 
two radio emission components in many sources in order to 
investigate the structure and evolution of the starbursts in 
these galaxies. 
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20 40 60 80 100 

t/Myr 

Figure 3. Top: Supernova rate as a function of time for a star 
formation rate that decays exponentially and a Salpeter IMF. 
The star formation rate has an e-folding time of t\> = 20 Myr 
and is shown as a dotted line. It is assumed that stars less mas- 
sive than 8 Mq do not produce type II supernovae; these have a 
lifetime ~ 50 Myr. The dashed lines show the supernova rate for 
stars still embedded in HII regions for various HII region lifetimes. 
From bottom to top these are, i HII = (10,20,30,40) Myr. Bot- 
tom: The fraction of the total supernova rate occurring within 
HII regions for the above values of tmi- We equate this fraction 
with the observed fractional radio emission of the more and less 
absorbed components, see table 2. 




HII region lifetime/Myr 

Figure 4. Fraction of supernovae occurring within dense HII 
regions as a function of the HII region lifetime, tmi- The star 
formation rate decays exponentially, t/>(t) oc e~ t / t b. The solid, 
dotted and dashed lines are for ages, t agG = (15, 30, 50) Myr. The 
black curves are for e-folding time, ij, = 15 Myr while the coloured 
lines are for = 30 Myr. 
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Figure 5. Effect of a rising star formation rate. The black curves 
are the same as those in fig. 4 while the red curves are for a star- 
burst that rises in time as ip(t) oc e t l th . The green lines show the 
effect of a flatter IMF for high mass stars (dn/dM oc Af -215 , 
Baldry & Glazebrook, 2003) for an exponentially decreasing 
burst. 

5.1 Relation between supernovae and ionized gas 

As we have seen, the form of the radio spectra shown in 
Fig. 1 imply that synchrotron emission commonly emerges 
from at least two distinct regions with quite different emis- 
sion measures. Such a situation could arise either because 
supernovae occur in two different environments or because 
cosmic rays are transported from one environment to the 
other after their acceleration in supernova shocks. Because 
higher energy electrons lose energy to synchrotron and in- 
verse Compton losses more rapidly than low energy elec- 
trons, the low energy electrons emit synchrotron photons for 
longer and can therefore be transported further from their 
source. If all supernovae occur in dense ionized media, con- 
sistent with regions of higher emission measure, the lower en- 
ergy electrons may be transported to regions of lower emis- 
sion measure while the more energetic particles end their 
synchrotron emitting lives still within the dense gas. 

If the transport of relativistic electrons within HII re- 
gions is well-described by a diffusion coefficient, D, then they 
can be transported a distance d ~ (Dt) 1 ^ 2 before ceasing 
to emit synchrotron radiation (Condon 1992). Although we 
do not know the value of D in our sources, it has been esti- 
mated for the starburst galaxy NGC 2146 by Lisenfeld et al. 
(1996) to be in the range 1 — 4 x 10 28 cm 2 s -1 . If relativistic 
electrons emitting at 244 GHz have a synchrotron lifetime, 
t s ~ 5 x 10 4 yr in the strong magnetic fields characteriz- 
ing these sources, then they could be transported a distance 
d ^ 80 pc. As giant HII regions have sizes of order 10 pc 
and relativistic electrons may be transported by methods 
more efficient than diffusion (outflows or winds) the escape 
of relativistic electrons from HII regions seems feasible. 

However, our radio spectra are characterized by quite 
distinct components, as opposed to a gradual changes in 
spectral index. Because the synchrotron lifetime is only a 
weak function of the emission frequency, r s oc u" 1 ^ 2 , the 
difference in lifetime of electrons emitting at 8.4 GHz and 
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at 244 MHz is only a factor of 6 for a given magnetic 
field strength. Given the inevitable inhomogeneity of these 
sources on parsec scales it seems extremely unlikely that 
distinct spectral breaks could be caused by the preferential 
escape of lower energy cosmic rays. 

It is far more likely, in fact expected, that supernovae 
occur in different environments. If the supernovae responsi- 
ble for the synchrotron emission are those with progenitors 
with masses ^ 8 M (Portinari, Chiosi & Bressan 1998) then 
synchrotron emission is generated for ~ 50 Myr (for solar 
metallicity) after a burst of star formation, ie. the lifetime 
of an 8 M star. Massive stars are born in dense molecular 
clouds, which they then ionize once luminous. Once these 
stars have formed, the parent cloud immediately stars to be 
disrupted, first by stellar winds and then by the first super- 
nova explosions, which will occur after ~ 3 million years for 
a 100 M star. 

The synchrotron emission component that shows the 
higher frequency turn-over, implying a higher emission mea- 
sure, should therefore originate from cosmic rays accelerated 
by the most massive stars, that end their lives while still em- 
bedded in their parent HII regions. The less absorbed syn- 
chrotron component has its origin in supernova explosions 
of less massive stars (down to 8 Mq) that out-live the parent 
HII region. 

In those sources with 2 emission components the life- 
time of HII regions is therefore bounded, on the low side, 
by the lifetime of the most massive stars, ~ 3 Myr, and on 
the high side by the lifetime of an 8 Mq star. Only if HII 
region lifetimes lie in this range will emission components 
with very different emission measures be seen. For sources 
in which the burst of star formation began more recently 
than 50 Myr ago the constraint on HII region lifetime be- 
comes, 3 Myr < tmi < tage- The observed fraction of emis- 
sion in each radio component constrains the lifetime within 
this range. Below, we attempt to estimate the mean HII 
region lifetime from the observed spectra. 

5.2 HII region lifetimes 

The HII region lifetime derived from the relative strength 
of the two radio emission components depends on various 
starburst parameters that are not known. These are: the 
IMF, the temporal evolution of the star formation rate and 
the synchrotron emission from a supernova as a function of 
progenitor mass. Bearing these limitations in mind, we can 
nonetheless make use of the modeling carried out in Vega 
et al. (2008). In this work, we modeled the spectral energy 
distribution from the near infrared to the radio of the present 
sample of LIRGs and ULIRGs with starburst models. For 
each object, both the starburst age, t agc , and the e-folding 
time of an assumed exponentially decaying burst, tb, were 
derived. 

Because the starburst timescales in these objects are 
thought to be very short, ie, of the the same order as the 
lifetime of the stars that undergo supernovae, we cannot 
assume that the supernova rate bears a one-to-one relation 
to the star formation rate. After the start of a burst, that 
then decreases in intensity with time, the supernova rate 
will initially rise. We estimate the supernova rate as follows. 
If the mass of a star and its lifetime are related via M(t), 
then for an IMF, tp(M), the number of stars as a function 



of their lifetime is ip(M(t)). If a starburst occurred at time 
t = 0, then at later time, t age , the number of stars about to 
end their lives is ip(M(t age —t)). The total supernova rate 
at time t age is then given by, 

/•tage-tmin 

SNR= cj>(t) i/>(M(t age - 1)) dt (2) 

Jti 

If t m in and t max are respectively the lifetimes of the most 
massive and least massive stars that end their lives as type 
II supernovae (t m in ~ 3 Myr for a 100 M star and t m ax ~ 
50 Myr for an 8 M star), then if t agc < t max then tl = 0, 
but if the burst is older than the lifetime of the longest 
lived supernova precursors, so that t agc > t max , then tl = 
(t agc — tmax). (j>{t) is the star formation rate. 

Some fraction of the supernovae, occurring at the time 
dependent rate given by equation 2, will explode while still 
embedded in their parent HII regions. These will be those 
stars with lifetimes shorter than the HII region lifetime, tmi- 
In equation 2 the lower limit on the integral becomes (t ag e — 
tmi) in this case. 

Because the synchrotron emitting lifetime of relativistic 
electrons, ~ 10 4 — 10 5 years in our frequency range, is much 
shorter than any timescale that characterizes the starbursts, 
we can interpret this supernova rate directly as the radio 
emission, if we assume that the radio emission per supernova 
is independent of the progenitor mass. 

We estimate M(t) from an analytical fit to the solar 
metallicity stellar evolution models of Bressan et al. (1993). 

log M = 1.389 - 0.351 log(t) + 0.223 log(t)~ 2 (3) 

We assume a Salpeter IMF, ip(M) oc Af~ 2 ' 35 , and ex- 
ponentially decreasing star formation rates, <f>(t) oc e~ t > lt *> , 
where tb is the burst e-folding time. Figure 3 illustrates how 
the total supernova rate varies as a function of time under 
these assumptions. The sharp peak is caused by the assump- 
tion that stars less massive than 8 M do not produce type 
II supernovae; these have a lifetime taken to be 50 Myr. 3 The 
dashed lines show the supernova rate for stars still embed- 
ded in HII regions for various values of tmi- We assume that 
all massive stars are born within HII regions and that after 
a time, tmi, these stars are no longer within their parent 
HII region. Stars with lifetimes shorter than tmi will ex- 
plode as supernovae while still within the HII region while 
those with lifetimes > tmi will produce supernovae external 
to HII regions. 



3 If stars with masses as low as, say, 5 M underwent type II su- 
pernovae the supernova rate would be unchanged up until 50 Myr 
after the start of a burst. Thereafter, the supernova rate would 
continue to increase until ~ 100 Myr, the lifetime of a 5 M star. 
The supernova rate peak would occur at 100 Myr rather than 
at 50 Myr. In both cases the fraction of embedded supernovae 
reaches a steady state after the peak. For a given HII region life- 
time, this steady value is a smaller fraction in the 5 M case 
(lower mass stars are more numerous and live longer, so the su- 
pernovae rate would be dominated by stars that have out-lived 
their HII regions). The effect of this, is that for sources with 
t ag c < 50 Myr, there is no change in the derived HII region life- 
times, whereas for those with t ago > 50 Myr the implied HII 
region lifetimes would increase. 
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By using the starburst e-folding time and starburst 
ages, as given by Vega et al. (2008), we can plot, for each 
source, the time variation of the total supernova rate and 
the fraction of supernovae occurring while still embedded 
in HII regions. The fraction of supernovae occurring in HII 
regions as a function of the HII region lifetime, tmi, for a 
given tage and t b is shown in Fig. 4. 

If we assume that the observed radio spectra reflect the 
properties of the mean HII region then these fractions can 
be compared with the observed radio flux ratio as given by 
the values of /2//1 in table 2. /1 then refers to radio emission 
from supernovae occurring within HII regions, while /2 refers 
to that produced outside HII regions. For each source, we 
therefore find the value of tmi, and these are given in table 4. 



6 DISCUSSION 

Individual HII regions cannot have lifetimes as long as most 
of those given in table 4. They can only remain ionized for 
a time equal to the lifetime of the least massive stars with 
sufficient ionizing flux. This is only a few Myr. Therefore, 
at least for an exponentially decreasing burst, the HII re- 
gions that cause the free-free absorption cannot be those 
that host most of the supernovae explosions that generate 
the radio emission. Younger star forming regions, still rich 
in ionized gas, may lie along the line of sight to a region 
of supernovae however, and cause the free-free absorption. 
Such a situation would naturally result if star formation were 
initiated in the nuclei of the mergers and then propagated 
outwards, perhaps triggered by supernova shocks. In addi- 
tion, this would create a shell-like (foreground screen) ge- 
ometry for the ionized gas and free-free absorption would 
be much more efficient. Therefore, what we actually con- 
strain is not the lifetime of an individual HII region in the 
usual sense, but rather, the time for which ionized gas, suf- 
ficient to cause the observed absorption, persists along the 
line-of-sight. This ionized gas need not be in a single 'HII 
region'. In what follows, we refer to 'regions of ionized gas' 
to avoid the connotation of a single cloud of gas. 

But how does the situation change if the starburst is 
in an increasing phase? For a rising star formation rate the 
required ionized gas lifetimes will decrease because the cur- 
rent star formation rate, resulting in embedded supernovae, 
is higher than that some Myr in the past, resulting in less 
embedded supernovae. Figure 5 compares the fraction of em- 
bedded supernovae as a function of tmi for an exponentially 
rising star formation rate, with one that falls exponentially, 
as in fig. 4. It is clear from this, that for a given fraction 
of embedded supernovae (absorbed radio flux), the required 
values of tmi are indeed much shorter. The rising star for- 
mation rate mitigates the accumulation of supernovae that 
occur outside regions of ionized gas, that otherwise cause 
the majority of supernovae to occur outside such regions. 
The effects of a rising star formation rate are important if 
t b < 50 Myr. 

Much shorter values of tmi would also result if the 
source ages were simply much lower than assumed. How- 
ever, given that the supernova rate does not peak until 
iage > 50 Myr, for almost any star formation rate and rea- 
sonable IMF, we cannot have a sample for which all the 
sources are younger than 50 Myr. We would have somehow 



selected out the older sources, that nonetheless would be 
strong radio emitters. Given that the original sample is ra- 
dio selected (Condon et al. 1991) the derived ages cannot be 
very much younger than those given in Vega et al. (2008). 

Another possibility is that the synchrotron emission is 
produced over shorter time scales. This could be achieved 
if the IMF were significantly flatter than the Salpeter value 
for high mass stars, as proposed by Baldry & Glazebrook 
(2003). However, the magnitude of the effect is very small, as 
shown in fig. 5. Shorter stellar lifetimes also result if source 
metallicities are higher than solar. For a metallicity 2.5 times 
the solar value (Z — 0.05) an 8 M Q star lives for 32 Myr, 
rather than 50 Myr, Fagotto et al.(1994). 

Until now we have assumed that the free-free absorption 
is caused by gas that has been ionized by the uv photons of 
young, massive stars. If this were not the case, then the time 
for which a region can remain ionized is no longer limited by 
the lifetime of massive stars. An AGN can, in principle, ion- 
ize large volumes of inter-stellar medium if it emits strongly 
at uv wavelengths and shorter. Vega et al. (2008) estimated 
the AGN contribution to the bolometric luminosity in the 
present sample of infrared luminous sources and found that 
only 9/30 sources have an AGN contribution that exceeds 
10%. Only IRAS 08572+3915 was found to be AGN dom- 
inated, with a 52% contribution from an AGN. Although 
AGN, where present, may contribute to the ionizing flux of 
a source, it seems unlikely that this would dominate over the 
contribution from star formation, which certainly produces 
ionizing radiation, and almost always dominates the AGN in 
terms of bolometric luminosity. An AGN that contributed 
additional ionizing flux, over and above that produced by 
star formation, would also create an excess of free-free emis- 
sion, visible as a flattening of the spectra towards higher 
frequencies. As mentioned above, and discussed in the next 
section, we actually find a deficit of free-free emission with 
respect to that expected for star formation alone. 

Finally, we note that the values that we derive for tmi 
are similar to the values of the 'escape time' found in Vega 
et al. (2008). This escape time parameterizes the time that 
elapses before stars escape the confines of the molecular 
clouds in which they were formed, and is described in Silva 
et al. (1998). The escape time is largely constrained by the 
spectral shape in the mid and far-infrared. Both these escape 
times and the values of tmi that we derive above, which are 
an order of magnitude larger than the values found in nor- 
mal star forming galaxies, suggest that the star formation in 
these objects is embedded in very large regions of dense gas. 
The effects of age-dependent extinction are thus also seen in 
the radio regime, as in other parts of the spectrum. 



6.1 Where is the thermal emission? 

The thermal fractions shown in table 2 are almost always 
lower than the 0.1 assumed by Condon (1992) for "nor- 
mal" galaxies. Despite this, there are still many instances 
where the data lie below the fitted curve. Among these, 
there are those that show a steepening of the spectra to- 
wards higher frequencies (eg. Arp 299), something that even 
a complete absence of thermal emission cannot achieve. The 
simple model we apply to these sources therefore fails at 
high frequencies. There is some additional effect influencing 
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the radio emission in these luminous infrared galaxies that 
is not required to fit the radio spectra of normal galaxies. 

The absorption of ionizing photons by dust can re- 
duce the quantity of ionized gas in ULIRGs (Valdes et 
al. 2005) thus reducing the thermal emission. However, al- 
though this effect may be present it cannot explain the spec- 
tra of those sources for which the spectra actually steepen 
above ~ 10 GHz. 

As discussed in Clemens et al. (2008), synchrotron ag- 
ing is not expected to cause a steepening of the radio spectra 
towards high frequencies in these sources because the syn- 
chrotron lifetime is very short. One actually expects spectra 
to be globally steeper than those of normal galaxies. The fact 
that this is not the case has been discussed by Thompson et 
al. (2006). These authors consider the effects of ionization 
losses on the radio spectra of starburst galaxies and find that 
the radio spectra gradually flatten towards low frequencies 
and that the spectral indices in the GHz range are similar 
to those of normal galaxies. If ionization losses are indeed 
important in these kind of galaxies it seems likely that the 
effects of ionization losses could explain why we observe un- 
expectedly steep spectra at high frequencies in fig. 1. 



7 CONCLUSIONS 

We have combined new measurements of the 244 and 
610 MHz fluxes, obtained with the GMRT, of a sample of 
20 LIRGs and ULIRGs with existing radio data. The result- 
ing radio spectra show a range of bends and turn-overs that 
can be well explained in terms of free-free absorption. In 
several cases two radio components with different emission 
measures are necessary to fit the data. 

For those objects with two emission components, we es- 
timate the fraction of radio emission associated with each. 
Via a simple model for the temporal evolution of the super- 
nova rate, we estimate the fraction of supernovae occurring 
inside and outside of their parent HII regions as a function 
of HII region lifetime. For simple exponentially decreasing 
star formation rates the lifetimes so derived are too long for 
them to refer to single HII regions. Altering the form of the 
high mass end of the IMF or the metallicity has relatively 
little affect on this result. However, if the star formation rate 
for these sources is currently rising, rather than declining, 
the derived HII region lifetimes are very much reduced. The 
effect of a rising star formation rate is to increase the num- 
ber of younger (high mass) supernovae occurring within HII 
regions, relative to the number of older (lower mass) super- 
novae occurring outside HII regions. The star formation rate 
rise times should be much less than the lifetime of the least 
massive stars that produce type II supernovae. Despite this 
reduction, the derived ages remain longer than those plau- 
sible for individual HII regions. Such large ages imply that 
the physical situation is one in which younger HII regions 
obscure older star formation regions along the line of sight. 
This situation would arise naturally if star formation prop- 
agated outwards during a burst. Because a rising star for- 
mation rate creates conditions more conducive to free-free 
absorption, and the fact that the number of sources with 
two emission components is close to 50%, it is likely that 
those in which only a single emission component is present 
are those with falling star formation rates. 



Table 4. Derived HII region survival times for a Salpeter IMF. T 
is the lesser of the source age as derived in Vega et al. (2008) and 
50 Myr; the lifetime of an 8 Mq star. tjj n are the survival times 
for exponentially decreasing star formation rates, for expo- 
nentially increasing star formation. In both cases the e-folding 
time for the star formation in each source is as derived in Vega 
et al. (2008). 



Source 


tage 


tesc 


T 


*HII 


t+ 

l HII 




Myr 


Myr 


Myr 


Myr 


Myr 


NGC34 


75 


65 


50 






1C1623 


43 


24 


43 


33 


25 


CGCG436-30 


32 


24 


32 






101364-1042 


19 


14 


19 






IIIZw35 


35 


35 


35 






UGC2369 


69 


36 


50 


>45 


>40 


103359+1523 












JNCj01o14 


69 


50 


50 


>42 


>30 


TAri on o c o a 


5.6 


6 


5.6 






NGC2623 


53 


30 


50 


43 


39 


108572+3915 


3.9 


3 


3.9 


1 


1 


T T/~1 f~i A OO 1 

UGC4881 


97 


44 


50 


>41 


>36 


UGC5101 


28 


22 


28 


23 


15 


I10173+0828 


16 


14 


16 


8 


8 


I10565+2448 


81 


60 


50 


>46 


>44 


Arpl48 


54 


27 


50 


>39 


>25 


UGC6436 


98 


34 


50 


>47 


>46 


Arp299 


57 


43 


50 


39 


27 


112112+0305 


7.6 


7 


7.6 


4 


4 


UGC8058 


37 


26 


37 


31 


27 


UGC8387 


52 


33 


50 


45 


39 


NGC5256 


43 


15 


43 






UGC8696 


45 


37 


45 


34 


22 


114348-1447 


15 


14 


15 


>12 


>11 


IZwl07 


59 


40 


50 


46 


44 


115250+3609 


19 


16 


19 






Arp220 


22 


20 


22 


16 


15 


NGC6286 


44 


23 


44 


>36 


>16 


NGC7469 


87 


45 


50 






IC5298 


79 


60 


50 






MRK331 


81 


35 


50 







The existence of two emission components in the radio 
therefore points to an age-dependent extinction analogous to 
that seen at other wavelengths. In the case of radio emission, 
the most massive stars, that end their lives earlier, produce 
synchrotron emission that is subjected to preferentially more 
free-free absorption than those, lower mass, stars that live 
for longer. 

The observed over-prediction of fluxes at frequencies 
above ~ 10 GHz may be resolved if the intrinsic synchrotron 
spectrum were not straight, but had a steeper spectral index 
at higher frequencies due to ionization losses, as modeled by 
Thompson et al. (2006). 

The two radio components in UGC 8058 can be identi- 
fied with the radio loud AGN and star formation; the former 
contributing 66% of the radio emission. This interpretation 
is consistent with the model of Vega et al. (2008). 
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